The effect of double (S238F/W159H) mutations on the structure and dynamics of PET degrading enzyme.

Polyethylene terephthalate (PET) is one of the highly produced synthetic polymers worldwide and had acquired attention due to its impact resistance, high clarity, and light weight. PET has become the first choice in making disposable bottles, leading to massive scales of production resulting in very high utilization across various facets of our daily life. Unfortunately, PET accumulates as waste and is highly resistant to biodegradation, thus presenting a serious threat to the ecosystem. Degradation of PET by enzymatic hydrolysis is a promising strategy to depolymerize the PET into its monomers. In recent studies, a plastic-degrading enzyme known as PETase (IsPETase) from the Ideonella sakaiensis has been identified to hydrolyze PET. The wild-type enzyme from Ideonella sp., has been engineered to improve the catalytic activity. While the IsPETase and its variants have been the subject of extensive structural and biochemical studies, the corresponding computational studies to support the improved activity of the mutant enzyme is not fully understood. In this work, we employed all-atom classical molecular dynamics simulations of the wild-type and double mutant IsPETase enzymes to investigate the underlying reason for the improved catalytic activity in the double mutant by means of structure-dynamics-function relationship. Our results show that the engineered mutations reshape the active site structure, volume, and dynamics of the protein loops which is crucial for substrate binding. We also demonstrate that addition of aromatic and hydrogen bond-forming residues near catalytic site improves binding affinity. This work will enable the rational design of mutants for enhanced PET degrading activity.Communicated by Ramaswamy H. Sarma.

Exploring the Conformational Impact of Glycine Receptor TM1-2 Mutations Through Coarse-Grained Analysis and Atomistic Simulations.

Pentameric ligand-gated ion channels (PLGICs) are a family of proteins that convert chemical signals into ion fluxes through cellular membranes. Their structures are highly conserved across all kingdoms from bacteria to eukaryotes. Beyond their classical roles in neurotransmission and neurological disorders, PLGICs have been recently related to cell proliferation and cancer. Here, we focus on the best characterized eukaryotic channel, the glycine receptor (GlyR), to investigate its mutational patterns in genomic-wide tumor screens and compare them with mutations linked to hyperekplexia (HPX), a Mendelian neuromotor disease that disrupts glycinergic currents. Our analysis highlights that cancer mutations significantly accumulate across TM1 and TM2, partially overlapping with HPX changes. Based on 3D-clustering, conservation, and phenotypic data, we select three mutations near the pore, expected to impact GlyR conformation, for further study by molecular dynamics (MD). Using principal components from experimental GlyR ensembles as framework, we explore the motions involved in transitions from the human closed and desensitized structures and how they are perturbed by mutations. Our MD simulations show that WT GlyR spontaneously explores opening and re-sensitization transitions that are significantly impaired by mutations, resulting in receptors with altered permeability and desensitization properties in agreement with HPX functional data.

EnzyDock: Protein-Ligand Docking of Multiple Reactive States along a Reaction Coordinate in Enzymes.

Enzymes play a pivotal role in all biological systems. These biomachines are the most effective catalysts known, dramatically enhancing the rate of reactions by more than 10 orders of magnitude relative to the uncatalyzed reactions in solution. Predicting the correct, mechanistically appropriate binding modes for substrate and product, as well as all reaction intermediates and transition states, along a reaction pathway is immensely challenging and remains an unsolved problem. In the present work, we developed an effective methodology for identifying probable binding modes of multiple ligand states along a reaction coordinate in an enzyme active site. The program is called EnzyDock and is a CHARMM-based multistate consensus docking program that includes a series of protocols to predict the chemically relevant orientation of substrate, reaction intermediates, transition states, product, and inhibitors. EnzyDock is based on simulated annealing molecular dynamics and Monte Carlo sampling and allows ligand, as well as enzyme side-chain and backbone flexibility. The program can employ many user-defined constraints and restraints and classical force field potentials, as well as a range of hybrid quantum mechanics-molecular mechanics potentials. Herein, we apply EnzyDock to several different kinds of problems. First, we study two terpene synthase reactions, namely bornyl diphosphate synthase and the bacterial diterpene synthase CotB2. Second, we use EnzyDock to predict reaction coordinate states in a pair of Diels-Alder reactions in the enzymes spirotetronate AbyU and LepI. Third, we study a couple of racemases: the cofactor-dependent serine racemase and the cofactor independent proline racemase. Finally, we study several cases of covalent docking involving the Michael addition reaction. For all systems we predict binding modes that are consistent with available experimental observations, as well as with theoretical modeling studies from the literature. EnzyDock provides a platform for generating mechanistic insight into enzyme reactions, useful and reliable starting points for in-depth multiscale modeling projects, and rational design of noncovalent and covalent enzyme inhibitors.

Chemical and Biochemical Approaches for the Synthesis of Substituted Dihydroxybutanones and Di- and Tri-Hydroxypentanones.

Polyhydroxylated compounds are building blocks for the synthesis of carbohydrates and other natural products. Their synthesis is mainly achieved by different synthetic versions of aldol-coupling reactions, catalyzed either by organocatalysts, enzymes, or metal-organic catalysts. We have investigated the formation of 1,4-substituted 2,3-dihydroxybutan-1-one derivatives from para- and meta-substituted phenylacetaldehydes by three distinctly different strategies. The first involved a direct aldol reaction with hydroxyacetone, dihydroxyacetone, or 2-hydroxyacetophenone, catalyzed by the cinchona derivative cinchonine. The second was reductive cross-coupling with methyl- or phenylglyoxal promoted by SmI2, resulting in either 5-substituted 3,4-dihydroxypentan-2-ones or 1,4 bis-phenyl-substituted butanones, respectively. Finally, in the third case, aldolase catalysis was employed for synthesis of the corresponding 1,3,4-trihydroxylated pentan-2-one derivatives. The organocatalytic route with cinchonine generated distereomerically enriched syn-products (de = 60-99%), with moderate enantiomeric excesses (ee = 43-56%) but did not produce aldols with either hydroxyacetone or dihydroxyacetone as donor ketones. The SmI2-promoted reductive cross-coupling generated product mixtures with diastereomeric and enantiomeric ratios close to unity. This route allowed for the production of both 1-methyl- and 1-phenyl-substituted 2,3-dihydroxybutanones at yields between 40-60%. Finally, the biocatalytic approach resulted in enantiopure syn-(3 R,4 S) 1,3,4-trihydroxypentan-2-ones.

Nucleoside-2',3'/3',5'-bis(thio)phosphate antioxidants are also capable of disassembly of amyloid beta42-Zn(ii)/Cu(ii) aggregates via Zn(ii)/Cu(ii)-chelation.

Currently, there is an urgent need for biocompatible metal-ion chelators capable of antioxidant activity and disassembly of amyloid beta (Aß)-aggregates as potential therapeutics for Alzheimer's disease (AD). We recently demonstrated the promising antioxidant activity of adenine/guanine 2',3' or 3',5'-bis(thio)phosphate analogues, 2'-dA/G3'5'PO/S and A2'3'PO/S, and their affinity to Zn(ii)-ions. These findings encouraged us to evaluate them as agents for the dissolution of Aß42-Zn(ii)/Cu(ii) aggregates. Specifically, we explored their ability to bind Cu(ii)/Zn(ii)-ions, the geometry and stoichiometry of these complexes, Cu(ii)/Zn(ii)-binding-sites and binding mode, and the ability of these analogues to dissolve Aß42-Zn(ii)/Cu(ii) aggregates, as well as their effect on the secondary structure of those aggregates. Finally, we identified the most promising agents for dissolution of Aß42-Zn(ii)/Cu(ii) aggregates. Specifically, we observed the formation of a 1 : 1 complex between 2'-dG3'5'PO and Cu(ii), involving O4 ligands. Zn(ii) was coordinated by both thiophosphate groups of 2'-dA3'5'PS and A2'3'PS involving O2S2 ligands in a 1 : 1 stoichiometry. A2'3'PS dissolves Aß42-Zn(ii) and Aß42-Cu(ii) aggregates as effectively as, and 2.5-fold more effectively than EDTA, respectively. Furthermore, 2'-dG3'5'PS and A2'3'PS reverted the Aß42-M(ii) structure, back to that of the free Aß42. Finally, cryo-TEM and TEM images confirmed the disassembly of Aß42 and Aß42-M(ii) aggregates by A2'3'PS. Hence, 2'-dG3'5'PS and A2'3'PS may serve as promising scaffolds for new AD therapeutics, acting as both effective antioxidants and agents for solubilization of Aß42-Cu(ii)/Zn(ii) aggregates.